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Abstract. In this talk, I briefly review several models of the heavy quarkonium production at 
collider energies, and discuss the status of QCD factorization for these production models. 

INTRODUCTION 

The production of bound states of heavy quark pairs has been the subject of a vast the- 
oretical literature and of intensive experimental study, and offers unique perspectives 
into the formation of QCD bound states [1]. The first step in quarkonium production, 
the inclusive creation of a pair of heavy quarks, is an essentially perturbative process, 
and takes place at a distance scale much smaller than the physical size of a quarkonium. 
Therefore, the transition from the produced heavy quark pairs to bound mesons is un- 
likely to be instantaneous. The "long" lifetime of the produced pairs allows the dynamics 
of QCD confinement to evolve, and provides us with a window of opportunities to probe 
the formation of QCD bound states via the interaction between the pairs and the medium 
where they were produced. 

Much of the predictive content of QCD perturbation theory is contained in factoriza- 
tion theorems [2]. They allow us to separate long-distance physics from short-distance 
interactions in hadronic cross sections, and to provide physical content for the uncalcu- 
lable long-distance quantities, so that they can be measured independently or calculated 
numerically. In this talk, I first briefly review several models of heavy quarkonium pro- 
duction, and then, discuss the status of QCD factorization for these production models, 
and finally, give a brief summary and outlook. 



PRODUCTION MODELS 

In order to produce a heavy quarkonium in hadronic collisions, the energy exchange 
in the collisions has to be larger than the invariant mass of the produced quark pair 
(> 2mg with heavy quark mass mg). The pairs should be produced at a distance scale 
Ar < l/2m,Q < 0.1 fm for charmonia or 0.025 fm for bottomonia. Since the binding 
energy of a heavy quarkonium of mass M is much less than heavy quark mass, (M 2 — 




4mg)/4mg <C 1, the transition from the pair to a meson is sensitive to soft physics. The 
quantum interference between the production of the heavy quark pairs and the transition 
process is powerly suppressed by the heavy quark mass, and the production rate for a 
heavy quarkonium state, H, up to corrections in powers of 1/mg, can be factorized as, 

OA+B^H+X ~ £ / dT Q Q A +B^QQ[n]+x( T QQ> m Q) F QQ[n]^H( T Qq) C 1 ) 

with a sum over possible QQ[n] states and an integration over available QQ phase 
space cITqq. Neglecting the "high twist" interaction between incoming hadrons and 
quarkonium formation, the nonperturbative transition probability F for a pair of off- 
shell heavy quark (y/) and antiquark (%) to a quarkonium state H is proportional to the 
Fourier transform of following matrix elements 

^\x\yi)^nW{y2)\H+N){H+NW\y 2 ).Kx{ym, (2) 

N 

where yi(yi) are coordinates, and J(f n and are local combinations of color and spin 
matrices for the QQ state n. A proper insertion of Wilson lines to make the operators 
in Eq. © gauge invariant is implicit [j]. hi Eq. ©, Ga+b-^qqM+x represents the 
production of a pair of on- shell heavy quarks and is calculable in perturbative QCD 
The debate on the production mechanism has been focusing on the transition from 
the pair to the meson. 

Color-singlet model The color-singlet model (CSM) assumes that only a color singlet 
heavy quark pair with the right quantum number can become a quarkonium of the same 
quantum number and the transition from the pair to a meson is given by the quarkonium 
wave function [5]. By neglecting the dependence on the pair's relative momentum in 
a A+B^QQ[n]+x^ I ^QQ^QQ^H m Eq. © is set to equal the matrix element in Eq. 
evaluated at y ; (y ; ) = 0, which is proportional to the square of coordinate-space quarko- 
nium wave function at the origin, \Rh(0)\ 2 [1], and therefore, 

OaIb-^h+x - G A+ B^QQ[H]+x(m Q ) \R H (0) | 2 • (3) 
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FIGURE 2. From Ref. @l, NRQCD predictions of charmonium polarizations are compared with the 
CDF data EH. 



The same wave function appears in both production and decay, and the model provides 
absolutely normalized predictions. It works well for J/i// production in deep inelastic 
scattering, photon production, and some low energy experiments yj], but fails to predict 
the CDF data, as demonstrated by the dotted lines in Fig.[T{a) @]. 

Color evaporation model The color evaporation model (CEM) takes a very different 
approach to handle the nonperturbative transition. It assumes that all QQ pairs with 
invariant mass less than the threshold of producing a pair of open-flavor heavy mesons, 
regardless their color, spin, and invariant mass, have the same probability to become a 
quarkonium L8]. That is, the FqqM-^h m Eq. © is a constant for a given quarkonium 
state, H, and therefore, 

pTMq 

0A+B-+H+X ~ /fl / dm QQ a A+B^QQ+x( m QQ) ( 4 ) 
J2mn 

with open-flavor heavy meson mass Mq and a constant /# Hfl. With a proper choice of 
m,Q and the model gives a reasonable description of almost all data including the 
CDF data, as shown in Fig.[TJb). 

Nonrelativistic QCD model The Nonrelativistic QCD (NRQCD) model is based on the 
fact that the typical heavy quark rest-frame kinetic energy and binding energy, mgv 2 , in 
a heavy quarkonium is much smaller than the heavy quark mass. The model separates 
the physics at scales of order mg and higher from the dynamics of the binding by using 
NRQCD, an effective field theory [ 10] . It provides a systematic prescription to calculate 
the physics at mg order-by-order in powers of a s , and expands the nonperturbative 
dynamics in terms of local matrix elements ordered in power series of heavy quark 
velocity v BaLLO], 

a A+B-?H+X = H°A+B^QQ[n]+X (&n) > ( 5 ) 
n 

where the are NRQCD operators for the state H ifToll. 

(0) = x\0) JT„ y(0) (a f H a„) ^(0) ,<z(0) , (6) 
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FIGURE 3. Sketch for process A+B — > H(p) +X: (a) sample scattering amplitude with possible 
factorization breaking soft interactions indicated by dashed lines, and (b) factorization with Wilson lines 
indicated by thin double lines. 

where a\ is the creation operator for H, % (i/a) are two component Dirac spinors, and 
and J^/ are defined in Eq. © and can also involve co variant derivatives. At higher 
orders in v, the operator 6^ can have additional dependence on field strength as well as 
more fermion fields. The factorization in Eq. © could be understood from Eq. dU) by 
expanding the GA+B->QQ[n]+x at heavy quark relative momentum, q = (pg — Pq)/2 = 0. 
The moments, J dT 'qq Q N Fqq^ h , lead to local matrix elements with high powers of v. 

The NRQCD model allows every QQ [n] state to become a bound quarkonium, while 
the probability is determined by corresponding nonperturbative matrix elements 
It in principle includes the physics of CSM. Its octet contribution is the most important 
one for high pj quarkonium production at collider energies [6]. The NRQCD model has 
been most successful in interpreting data |H0], as shown in Fig.[TJa). 

Quarkonium polarization and other models The key difference between the NRQCD 
model and the CEM is the prediction on quarkonium polarization. Once the matrix el- 
ements are determined, the NRQCD model can systematically calculate the po- 
larization of produced heavy quarkonia. On the other hand, CEM does not provide a 
systematic prescription to calculate the polarization because of the uncontrolled radi- 
ations from the quark pair. The polarization of quarkonia at large pj was considered 
a definite test of the NRQCD model [6]. But, the model has failed the test, as seen in 
Fig. El if the CDF data hold H 1 ill - Several improved and new models have been proposed 
to address the issues of polarization ifl^l . 



FACTORIZATION 

The heavy quarkonium production in hadronic collisions involves both perturbative and 
nonperturbative scales. The nonperturbative physics appears not only in the transition 
from the heavy quark pair to a bound state but also in incoming hadron wave functions. 
A typical scattering amplitude for quarkonium production, as sketched in Fig. Eta), 
can have soft and nonperturbative interactions between incoming hadrons as well as 
between the spectators and the formation process. These soft interactions may introduce 
process dependence to the nonperturbative matrix elements, and consequently, spoil the 



predictive powers of Eqs. @ and ©. 

A proof of the factorization needs to: 1) show that the square of the scattering 
amplitude in Fig. Eta), after summing over all amplitudes with the same initial and final 
states, can be expressed as a convolution of the probabilities, as sketched in Fig. 0b); 
each probability represents a square of sub-amplitudes and is evaluated at its own 
momentum scale(s); 2) show that the piece evaluated at perturbative scale(s) is infrared 
safe and those evaluated at nonperturbative scales are universal. 

As argued in Ref. 0], cross sections for producing on-shell heavy quark pairs can be 
computed in terms of QCD factorization. Therefore, the right-hand-side of the CEM 
formalism in Eq. © can be factorized into a convolution of a calculable partonic 
cross section of producing the heavy quark pair and two universal parton distributions 
from respective hadrons. However, the factorization statement here does not provide 
justification that the FqqM-^h m Eq. © is independent of the pair's invariant mass 
kiqq, spin, and other quantum numbers. For quarkonium production at a large pr, the 
should be universal within the model because soft interactions in Fig.|3ta) are suppressed 
by powers of l/pr- When pj <C niQ, the fn may not be universal. 

Fully convincing arguments have not yet been given for NRQCD factorization for- 
malism in Eq. © BfALLO]- Since the interaction between the beam jet and the jet of heavy 
quark pair should be suppressed by powers of l/pr, one might expect the NRQCD 
factorization formalism to work at large pj- When pj ^> mg, the heavy quarkonium 
production is similar to the single light hadron production, and is dominated by par- 
ton fragmentation. The cross section is proportional to the universal parton-to-hadron 
fragmentation functions yD, 

A+b ^h+x(pt) = £ &A+B-i+x(PT/z,ll) ®D H/i (z,m Q ,n) + @{m 2 H /p\) . (7) 

Here, ® represents a convolution in the momentum fraction z- The cross section 
&A+B-*i+x includes all information on the incoming state, including convolutions with 
parton distributions of hadrons A and B at factorization scale /i, as sketched in Fig. 01b). 
As a necessary condition for NRQCD factorization in Eq. ©, the following factorization 
relation, 

D H/i (z,m Q ,n) = £ ^ 2 2 [n] (z,jU,«e) (O , (8) 

n 

is required to be valid to all orders in a s and all powers in v-expansion for all parton- 
to-quarkonium fragmentation functions Jj]. In Eq. ©, ^i-^QQ[ n ] describes the evolution 
of an off-shell parton into a heavy quark pair in state [«], including logarithms of ji /niQ, 
and should be infrared safe [3|]. 

The factorization relation in Eq. (0) was tested up to next-to-next-to-leading order 
(NNLO) in a s at v 2 order in Ref. 0], as well as at finite v in Ref. lfl3ll . Consider 
representative NNLO contributions to the fragmentation process of transforming a color 
octet heavy quark pair to a singlet, as sketched in Fig. |4j The individual classes of 
diagrams in Fig. |4jl) and (II), for which two gluons are exchanged between the quarks 
and the Wilson line, satisfy the infrared cancellation conjecture of Ref. [10], by summing 
over the possible cuts and connections to quark and antiquark lines, as do diagrams that 
have three gluon-eikonal vertices on the quark pair and one on the Wilson line Jj]. For 
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FIGURE 4. Representative NNLO contributions to g — > QQ fragmentation in eikonal approximation, 
see Ref. |3] for the details. 



Fig. EtlH) type of diagrams, however, with a three-gluon interaction, this cancellation 
fails. By summing over all contributions to second order in the relative momentum q, it 
was found that the fragmentation function has a noncanceling real pole fl] 



2 1 2 

s 3e 



(9) 



From Fig. Eflll), this infrared divergence is not topologically factorizable and can not 
be absorbed into the matrix elements As demonstrated in Ref. JH], as a necessary 

condition for restoring the NRQCD factorization at NNLO, the conventional & (v 2 ) octet 
NRQCD production matrix elements (&„) must be modified by incorporating Wilson 
lines that make them manifestly gauge invariant, so that the infrared divergence in Eq. © 
can be absorbed into the gauge-completed matrix elements. 

For quarkonium production, all heavy quark pairs with invariant mass less than 
a pair of open-flavor heavy mesons could become a bound quarkonium. Therefore, 
the effective velocity of heavy quark pair in the production could be much larger 
than that in decay. For charmonium production, charm quark velocity, v c ~ \q c \/ m c < 



(4Mp — 4m 2 ) / (4m 2 ) ~ 0.88, is not small, and therefore, the velocity expansion for 
charmonium production may not be a good approximation, unless one can identify and 
resum large contributions to all order in v or have a factorized formalism at finite v. 
It was found in Ref. lfl3ll that with the gauge-completed matrix elements, infrared sin- 
gularities in the fragmentation function for a color octet pair to a singlet at NNLO are 
consistent with NRQCD factorization to all orders in v or to a finite v, 
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(10) 



where f(v) = 2v/(l + v 2 ). The result in finite v is remarkably compact and intriguing, 
and should encourage further work on the factorization theorem. 



SUMMARY AND OUTLOOK 



After more than 30 years since the discovery of J/i//, we still have not been able to fully 
understand the production mechanism of heavy quarkonium in high energy collisions, 



in particular, the transition from the heavy quark pair to a bound quarkonium. Although 
there are good reasons for each production model, none of the factorized production 
formalism, including that of the NRQCD model, has been proved theoretically. Further 
work on the factorization theorem is critical. 

The heavy quarkonium production has a "long" transition time from the produced 
pairs to bound mesons, and can be a good probe for the properties of newly formed hot 
and dense matter of quarks and gluons, or the quark-gluon plasma, in ultrarelativistic 
heavy ion collisions ffil4L ilSfl - if we can calibrate the production. Since the medium in- 
teraction is sensitive to the quantum state of the proporgating heavy quark pair, nuclear 
matter could be an effective filter to distinguish the production mechanism [16]. De- 
tailed studies of nuclear dependence of heavy quarkonium production in hadron-nucleus 
collisions should provide invaluable information on the formation of heavy quarkonia in 
hadronic collisions. 
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